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Abstract 

The renormalization group for maximally anisotropic su(2) current interactions in 2d is shown to 
be cyclic at one loop. The fermionized version of the model exhibits spin-charge separation of the 
4-fermion interactions and has symmetry. It is proposed that the S-matrices for these theories 
are the elliptic S-matrices of Zamolodchikov and Mussardo-Penati. The S-matrix parameters are 
related to lagrangian parameters by matching the period of the renormalization group. All models 
exhibit two characteristic signatures of an RG limit cycle: periodicity of the S-matrix as a function 
of energy and the existence of an infinite number of resonance poles satisfying Russian doll scaling. 
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I. INTRODUCTION 



Any solution of the Yang-Baxter equation that is crossing-symmetric and unitary be- 
comes a candidate for the factorizable S-matrix of a quantum field theory in 2 space-time 
dimensions Q|. Though infinite classes of rational and trigonometric solutions related to 
quantum algebras are known to describe an underlying quantum field theory, there are no 
a priori principles that ensure that every S-matrix necessarily has such an interpretation. 
In this regard the elliptic solutions provide a striking example; though these solutions have 
been known for over 25 years, it has remained unknown whether they have a field the- 
ory description. In this paper a field theory will be proposed for Zamolodchikov's elliptic 
S-matrixj^j], and also for the Mussardo and Penati's scalar S-matrix^]. 

An elliptic solution to the Yang-Baxter equation was first found by Baxter in connection 
with the 8 vertex model of classical lattice statistical mechanics (^j). It was later shown that 
it could be made crossing symmetric and unitary 2]. We emphasize that the answer to the 
question of whether there is a field theory description underlying this S-matrix is not con- 
tained in the physics of the 8 vertex model, nor equivalently the fully anisotropic XYZ spin 
chain. There is no a priori relation between the solution of the Yang-Baxter equation that 
encodes the Boltzman weights and the S-matrix of the fundamental excitations. Further- 
more the XYZ spin chain has a non-relativistic dispersion relation, whereas the S-matrix is 

relativistic. Certain relativistic continuum limits of the XYZ spin chain exist that lead to 

f 

the sine-Gordon, or massive Thirring model [5], but the resulting S-matrix is not elliptic but 
trigonometric. 

The new ingredient we use to relate the elliptic S-matrices to a field theory is the limit- 
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cycle in the renormalization group (RG) flow of the field theory 

Indeed, it was suggested by Zamolodchikov|2( that the field theory, if it exists, should be 
characterized by an RG limit-cycle. His argument was based on the periodic properties of 
the S-matrix as a function of energy (at high energy) and the early arguments of Wilson 
showing that such periodicities are naturally a consequence of a cyclic RGy . (See section 
II.) 

An outline of this article is as follows. In section II we describe some physical signatures 



of an RG limit cycle in a general, model-independent way, extending the discussion in 
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to include the role of resonance poles. An infinite sequence of resonance poles leading to 
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masses with special scaling relations is another characteristic signature of an RG limit cycle, 
termed "Russian doll" scaling in [jj. There is a clear distinction between massive and 
massless theories. We argue that a massive theory can only correspond to a limit cycle in 
the ultra-violet, whereas a massless theory can support a limit cycle in the infra-red and 
ultra-violet, in fact on all scales. In section III we define the quantum field theory models 
as fully anisotropic su(2) current- current interactions in 2D. The bosonized form of the 
action involves both a scalar field and its dual. The fermionized version continues to exhibit 
spin-charge separation and has an explicit TL\ symmetry. 

In section IV we show that integration of the 1-loop beta functions leads to couplings 
expressed in terms of elliptic functions, and are hence periodic as a function of scale with 
the period a function of RG invariant combinations of the couplings. By analyzing the U(l) 
invariant limit we classify 4 possible models. In section V the invariant elliptic solutions 
of the Yang-Baxter equation are reviewed. These elliptic solutions are related to the field 
theory in sections VI and VII. In the U(l) invariant limit these theories go over to the usual 

n Qn 

sine-Gordon model[l|] or to the cyclic sine-Gordon model studied in L 9]|.37 ( j. That the elliptic 
S-matrices can have two different U(l) invariant limits relies on the fact that they have 
two different trigonometric limits at complementary elliptic moduli. The relation between 
lagrangian and S-matrix parameters is found by matching the period of the RG with the 
periodic S-matrix properties. The models also all have the expected Russian doll spectrum 
of resonances. There is another regime that is an elliptic deformation of the sinh-Gordon 
model and this is described in section VIII, where it is proposed that the S-matrix is the 

n 

scalar one considered by Mussardo-Penati|3|. Finally in section IX we consider the massless 
theories. In all cases the U(l) invariant limit serves as a non-trivial check. 

II. PHYSICAL SIGNATURES OF RG LIMIT CYCLES 

In this section we describe in a general, model-independent fashion some of the signatures 
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of an RG limit cycle, extending the discussion in 

Before proceeding, we mention that naively, limit cycle behavior in a unitary theory 
appears to be inconsistent with Zamolodchikov's c-theorem: c is a function of the running 
coupling constants, so if the couplings are periodic as a function of scale, so is c. This 
leaves us with a puzzle since our model is a hermitian perturbation of a unitary conformal 
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field theory. The S-matrix is also unitary in a strict sense. This issue was discussed at 
length in |l2| for a limiting case of the models here. Though further investigation is needed 
regarding this issue, we offer the following hint as to what may be going wrong. Instead of 
Zamolodchikov's c-function which is related to a two-point function of stress-energy tensors, 
let us consider the c-function c e g that determines the finite size effects; this is the function 
that was studied in For the quantum field on a cylinder of circumference R } the one- 
point function of the trace of the stress-energy tensor is a derivative of c e g: 

W* = ™ « 

It trivially follows that if T£ is positive, then c e g decreases with increasing R. As we now 
argue, the positivity of can generally be violated in theories with no ultra-violet fixed 
point and marginal perturbations, which is precisely the situation for our models. Suppose 
a quantum field theory is described by a conformally invariant action S c { t perturbed by 
operators O a \ 

S = S cit + J2 [p9AO A (2) 

where qa are positive couplings. Let (3a be the beta function for qa and Ya the scaling 
dimension of O a . Then the beta functions (for increasing length scale) to lowest order are 

p A = (2-T A )gA + 0(g 2 ) (3) 

The main point is that the trace of the stress-energy tensor recieves quantum corrections, 
and since it must be zero at a fixed point where the beta functions are zero, one must have: 

t^ = J2M 9 )O a (4) 

A 

The above formula is well known [22] and is easily verified to lowest order in conformal 
perturbation theory. Consider first a theory that can be formulated as a perturbation of 
an ultra-violet fixed point by relevant operators, which implies Ta < 2. Then the beta- 
functions are positive to lowest order. Furthermore, for relevant perturbations, because of 
the anomalous dimensions of the couplings, there is often no higher order corrections to 
the beta functions since higher powers of $a do not have the right dimension. So in this 
situation, is generally positive and the c-theorem holds. 

The above arguments clearly point to the way in which the c-theorem can be violated. 
If the O a are marginal, Y a = 2, and the beta functions start at 0(g 2 ). There are no 



general constraints on the sign of the beta function. For instance, suppose the theory has 
no ultra-violet fixed point so that the couplings increase toward the ultra-violet rather than 
go to zero. The latter implies the beta function must be negative. This in turn implies 
the stress-energy tensor need not be positive and the c-theorem is violated. For the models 
below, one can check that the beta functions are not always positive. 

Let us return to the properties of theories with RG limit cycles. Generally speaking, 
cyclic behavior in the RG flow can in principle exist at all scales, or can be approached 
asymptotically in the ultra-violet (UV) or infra-red (IR), the latter being UV or IR limit- 
cycle behavior. Once the flow is in the cyclic regime, the couplings are periodic 

g(l + \) = g(l) (5) 

where I = logL and L is the length scale. Above, the period A is fixed and model-dependent. 
It is well-known that the RG leads to scaling relations for the correlation functions, 



251 ]. Let G be an n-point correlation 



generally referred to as Callan-Symanzik equations 
function: 

G(xx,x 2 , .., x n ; g, I) = • • • $(ar n )) (6) 

The above correlation function is computed using the renormalized action and is thus finite 
and depends on the RG scale I. We have not explicitly displayed dependence on a mass 
parameter since for our models the action contains no such parameter: the physical mass m 
of particles is generated dynamically and is a function of g, I. The Callan-Symanzik equation 
expresses the independence of G on the arbitrary scale I: 

{Wl + m ifg + ni{9) ) G{X] 9J) = ° (7) 
where "f(g) is the anomalous dimension of $. Let dq> denote the naive (engineering) dimension 
of Then since G has dimension nd$>, the above equation leads to the scaling equation: 

'Ts + ^ 9) dg~ + n{l{9) ~ ^ ) ° {eSX] 9,l) = ° (8) 
The above equation can be explicitly integrated: 

G(e s x;g,l) = e~ snd * exp Uj^dg G(x;g(s),l) (9) 

where g(s) flows according to the beta-function: dg/ds = (3(g) with g(0) = g. Letting now 
s = A, one finds 

G(e x x;g,l) = e- Xnd *G(x;g,l) (10) 



Therefore the function G(x)/\x\ nd * is a periodic function of log |x| with period A. 

The S-matrix is related to the Fourier transform of the Green's functions to momentum 
space, and thus also obeys scaling relations. Let S denote the 2-particle to 2-particle S- 
matrix. For an integrable quantum field theory in 2 space-time dimensions, S only depends 
on the kinematic variable E 2 m = {Pi + P 2 ) 2 where Px,2 are the energy-momentum vectors 
of the incoming particles (for an integrable theory, the incoming and outgoing momenta are 
the same). Since the S-matrix is a dimensionless quantity one obtains: 

S(e- s E cm ,g) = S(E cm ,g(s)) (11) 

If a theory is in a cyclic regime, when s = A the above equation implies a periodicity in 
energy: 

S(e- x E cm ,g) = S(E cm ,g) (12) 



A. Ultra-violet limit cycles 

Let us now consider a UV limit cycle. We specialize to 2d kinematics, and first assume 
the spectrum of particles is massive, such that the energy-momentum can be parameterized 
in terms of a rapidity /3: 

E = mcosh./3, p = msmhf3 (13) 

Above, m is the physical mass of the particles, and as explained above, is an (unknown) 
function of g, I. The center of mass energy is 

E 2 m = 2m 2 (l + cosh/3), (3 = Pi — P2 (massive case) (14) 

The UV limit corresponds to high energies where j3 is large and E cm « me 13 / 2 . The relation 
eq. (|T2"j) then implies a periodicity in rapidity: 

S(P-2\) = S(0) (15) 

The above periodicity is the primary signature of a UV limit cycle for a massive theory. 

There can exist another signature in the UV that has analogies with properties of other 
models with IR limit cycles US, Namely, if {E n , g, L} is the spectrum of eigenvalues 

of the hamiltonian for a system of size L, then 

{E n ,g(s),e s L)} = {E n ,g,L}, (16) 



When s equals the period of one cycle A, the above equation shows that the energy spectrum 
at fixed g should reveal a discrete self-similarity as a function of log L. The manner in which 
the spectrum can reproduce itself after one RG cycle is dependent on the existence of an 
infinite number of eigenstates with the "Russian doll" scaling behavior: 



e x E n , n = 0,1,. .00 for n large (UV) (17) 



In each cycle these eigenstates reshuffle themselves such that the (n + l)'th state plays the 
same role as the n'th state of the previous cycle. Since we are considering a UV limit cycle, 
the above property holds at high energies E n . 

For an integrable quantum field theory in infinite volume, the above property can be 
manifested as an infinite tower of resonances with masses that obey special scaling relations. 
Namely, one expects resonances of mass M n , n = 0, 1, ..00, where 

M n+ i w e A M n , for n large (18) 

These resonances correspond to poles in the S-matrix. The physical strip is the region 
< Im{pi) < 7r, and the only allowed poles in this region correspond to stable bound states 
and are required to be on the imaginary axis. In the sequel, this requirement will provide 
some constraints on the models. Poles on the strip — 71 < Jm(/3) < with a non-zero real 
part correspond to unstable resonances. Other non-cyclic models with a finite number of 
lied in S-matrices with an infinite number of resonances 



resonances were studied 

were considered in 0, 0, Q . Consider a pole at (3 = fi — irj with \i > and < 77 < tc in 

the S-matrix for the scattering of two particles of mass m. This corresponds to a resonance 

of mass M and inverse lifetime V where 

2 

2 



M- — \ = 2m 2 (l + cosh(/i-ir/)). (19) 



Equivalent ly: 



r 2 

M 2 — — = 2m 2 (1 + cosh/i cos 77), 



Mr = 2m 2 sinh/i smr/. (20) 

In order for M, T to both be positive, both /i and 77 must be positive. Consider an infinite 
sequence of resonance poles at (3 n = /i n — irj n . When /i n is large one finds: 

M n w me' 1 " /2 cos(r/ n /2), T n w 2me^ /2 sin(^ n /2) (21) 



Comparing with eq. (JTBJ) one finds 

fi n+1 -fi n = 2X (UV) (22) 

The above equation relating the location of the real parts of the poles of the S-matrix is 
another primary signature of a UV limit cycle. 

Note that the above two signatures of a UV limit cycle are not independent. If the S- 
matrix has resonance poles and is also periodic in rapidity, then the 2A periodicity in rapidity 
automatically implies that the resonance poles satisfy eq. (122(1 . Our elliptic S-matrices will 
indeed have both the periodicity in rapidity and the Russian doll spectrum of resonances. 



B. Infra-red limit cycles 

Suppose the theory has a limit cycle that is approached in the infra-red rather than the 
ultra-violet. For now let us continue to suppose the spectrum is massive, though as we will 
argue, this is does not appear to be consistent. Here the periodicity in energy eq. (JT2J), does 
not directly lead to a periodicity in (3 since f3 is small at low energies and E cm ps 2m. Thus, 
for a massive theory, periodicity in rapidity of the S-matrix is not a signature of an IR limit 
cycle. 

Another possible signature is again based on eq. (|T5|) with s = X, which leads to a 
Russian doll scaling spectrum at low energies: 

E n+1 ps e~ x E n , for n large (23) 

Note the minus sign in comparison to the UV signature eq. ([17)1 . This corresponds to an 
accumulation of resonances near zero energy. Indeed the extension of the BCS hamiltonian 
in Q and the m „del in B has an IR llm , c y cle with the property e, 0. It turns out it 
is not possible to obtain a spectrum of masses scaling like in eq. fl23j) from resonance poles 
in the S-matrix as eq. (J20|) shows: cosh/i behaves always as e M with // > 0. 

Based on the above discussion we conclude that a massive theory cannot support a limit 
cycle in the IR. This is in accordance with the effective central charge computations in 
which show that c e ff is only quasi periodic in the UV and decays to zero in the IR for a 
massive theory. 
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Both the above problems are fixed if the theory is massless. Here the massless dispersion 
relations E = ±p can be parameterized as: 

7-1 m a 171 a 

E = —e PR , p = —e PR 1 tor right — movers 

TTt TTt 

E = — e _/3i , p— e _/3i , for left — movers (24) 

where now m is an energy scale. The center of mass energy for a right-mover with rapidity 
(3r scattering with a left-mover of rapidity (3l is 

Ecm = me 13 / 2 , (3 = (3r — /3l, (massless case) (25) 

If Srl(P) is the S-matrix for the scattering of right-movers with left-movers, then eq. (|12|) 
again implies a periodicity: 

S RL ((3-2X) = S RL (P) (26) 

One sees then that the main difference between the massive and massless case is that 
in the massive case the periodicity in rapidity eq. (|15|) implies a periodicity in energy E cm 
only for large E cm , whereas in the massless case it leads to a periodicity at all energy scales 
because E cm oc e 13 ^ 2 . Thus a massless theory with the periodicity eq. (|26J) is consistent with 
a cyclic RG flow on all scales, i.e. both the IR and UV. 

The same conclusion is reached when one considers resonances. Consider a pole in Slr{(3) 
at /3 = [i — it). In the massless case eq. (jTT)j) becomes: 

r\ 2 

M — i— J = m 2 exp(/i — irj) (27) 

and equation (|21j) is exact. The point is that unlike the massive case, now \l is allowed to 
be negative. An IR spectrum of masses satisfying 

M n+1 = e~ x M nj n = 0,l,...,oo, (IR) (28) 

is possible with an infinite sequence of resonance poles satisfying: 

- = -2A, {IR) (29) 

A massless model with resonance poles of real part 2Xn for all n positive or negative has 
both the IR and UV signatures of a limit cycle and can thus describe a cyclic RG on all 
scales. The massless S-matrices we consider in the sequel have this property. (See section 
IX.) 



III. THE MODELS AND THEIR BOSONIC AND FERMIONIC DESCRIPTIONS 



A. Fully anisotropic current perturbations 

We consider a conformal field theory with su(2) level k current algebra svmmetrv|l4l.lli 
perturbed by fully anisotropic current interactions: 

s = s clt + J^-a ( 9x rT + gyjyy + 9z j z T) (30) 

where g x 7^ g y 7^ g z are marginal couplings. Above, S c & is formally the action for the 
conformal field theory, e.g. the Wess-Zumino-Witten action at the critical point. The 
currents are normalized to have the following OPE: 

J a (z) J b (0) ~ ^5 ab + -J abc J c (0) (31) 

where z = (t + ix)/y2,z — (t — ix)/y/2 are euclidean light-cone space-time variables. 
(Minkowski space with real time is obtained by t — > it.) The structure constants are f abc = 
ie abc , where e abc is the completely antisymmetric tensor, and similarly for the right-moving 
currents J a {z). 

Our subsequent analysis will involve consideration of the U(l) current with components 
J z , J Z . We normalize this current as follows: 

j„ = (j z ,h) = ^(J,J) (32) 

and the U(l) charge T as 

T = J dx h = J dx (jz + h) (33) 

With this normalization, the currents J 1 * 1 = (J x ± iJ y )/\/2 have U(l) charge ±1. 

When the level k — 1, there are simple realizations of the current algebra in terms of a 
free boson or a doublet of free fermions 3] , which are described in the next subsections. 



B. Bosonization 

When k — 1 the su{2) currents have a free massless boson representation with Virasoro 
central charge c = 1. The action S c n is just the massless Klein-Gordon action. We normalize 
this action as follows: 

10 



so that the propagator is: 

(<j>(z,z)<j)(p)) = -logzz (35) 

In the conformal limit the free boson <fi can be separated into its right and left moving parts, 
<p = (f(z) + fp(z), and the currents have the following expressions: 

1 i 
J± = - 7 =exp(±zv / 2^), J z = ^=d z (p 

J = — exp(TiV2ip), f = -—djp (36) 
In this bosonic representation, the U(l) current is topological. In Minkowski space: 

where e 01 = — e 10 = 1, and it is identically conserved: d^j^ = 0. 
One finally finds for the action: 

S = ~L J d2x \{ {l + + 8g + cos ^ + 8g - cos ^) ( 38 ) 

where the dual field <f) = <p — Tp, and g± = g x ± g y . Noting that d z (p = d^p, &z4> = — <3z0, the 
relation between and its dual can be expressed in Minkowski space as: 

= e^d u (f) = 2V2n f (39) 

The operators exp(±zy / 20) have U(l) charge ±2 and thus the U(l) symmetry is broken 
when g x ^ g y . When g x = g y the dual field does not appear, and the model is the sine- 
Gordon model with T the usual topological charge for the U(l) symmetry Models similar to 



the one defined in eq.(38) have been studied in references QQI^Q) 21|, especially the self 
dual cases where g + = g- and when the perturbations of the gaussian model are relevant. 
In our model the latter perturbations are marginal and, as we shall see in section IV, the 
self dual constraint g + = g_ (e.g. g y = 0) is preserved by the RG only if g x or g z vanish. 
Both cases correspond to well known gaussian models. 



C. Fermionization 



For potential applications to condensed matter physics, we now consider a fermionic 
representation of the model. Again when k — 1, the su(2) currents can be represented as 

II 



fermion bilinears with a doublet of fermions in the spin 1/2 representation of su(2). Unlike 
the bosonic representation above, this is not an irreducible representation of the affine Lie 
algebra su(2) and instead has Virasoro central charge c = 2. It is thus a different model than 
that of the previous subsection. However we will show that they share the same S-matrices 
in the interacting sector. 

Introduce su(2) spin-1/2 doublets of left and right moving fermions: 

-(£)• *«-(::) 

and their hermitian conjugates where e.g. & L = (ip^, V'ij.)- The currents have the 

following representation: 

J± = \^ L a^ L , J z = \^ L a z * L (41) 

and similarly for J with L — > R, where <7j are the standard Pauli matrices and a ± = 
(cr x ± ia y )/2. The conformal action is now: 

S <* = f ^ E tylaMLa + lhd*1>Ra) (42) 

J o=T,4. 

The interaction terms are 

/d 2 x r 

The U{\) current is now: 

T= I^J dX (^T^T - ^L^Li + ~ ^Ri^Rl) (44) 

With this normalization the fields have the following charges: 

T = +1 : ipLi^Ri^lpipRi 

T = -l: fai^m^H^m (45) 

As before, using this one sees that the (g x — g y ) terms have charge T = ±2 and break 
the U(l) symmetry. However there is a remaining Z4 symmetry. Let % denote a finite 
U(l) transformation by an angle 9. Namely, if an operator O q has U(l) charge q, then 
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%{O q ) = e iq9 O q . Then, based on the charges in eq. (|45|h one finds %(J ± J 



,±4i0 



Thus the action is only invariant for 9 = 7r/2, and since (T^) 4 = 1, this corresponds to a 
ZU symmetry. The elliptic S-matrices in 3], which we will propose in the sequel to describe 
the model, where shown there to be Z4 symmetric. 

The fully anisotropic model continues to enjoy the so-called spin-charge separation. (For 
a discussion in the £7(1) invariant case, see for instance Q]) The fermions can be bosonized 
with two bosons d>t,d>\: 



(46) 



and the same with Defining bosons for the spin and charge degrees of freedom: 



V2 



V2 



(0T + 01 ) 



then the action is: 
S = — I ' d 2 x 

47T 



-(1 + 4g z )(d<j) s ) 2 + ^(<90 c ) 2 + Ag + cos yftfa + 4#_ cos yftfa 



(47) 



(48) 



Thus the spin and charge fields are decoupled. The C field is a free boson and the 4> s 
field has the same action as the boson of the c = 1 model eq. (|H8j) . Thus the S-matrices 
we propose below for the bosonic model also describe the scattering of the spin degrees of 
freedom for the fermionic model. This explains the Z4 symmetry of this S-matrix, though 
this symmetry is hidden in the bosonic description. 

IV. RENORMALIZATION GROUP AND CLASSIFICATION OF PHASES 
A. Beta functions and periods 



More generally, consider a conformal field theory perturbed by marginal operators O a : 



(49) 



Assuming the perturbing operators form a closed operator product expansion in the confor- 
mal theory: 

(50) 



A (z,z)0" (0,0) ~ —Cf? O°(0,0), 

'Z/'Z/ 
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the one-loop beta function is known to depend only on the coefficients C[22j: 

fo = %r = -\Y, C F9*9c (51) 

B,C 

where I, the RG 'time', is the log of the length scale, and the flow toward the infra-red 
corresponds to increasing I. For the current-current interactions defined in eq. (|3U|) . using 
eq. (|31|) one finds that the only non-zero C's are 

C xy = C yx = C zx = Q xz = Q yz = gzy = _ 4 ^ 

This gives 

13 x = 4g y g z , (3 y = Ag x g z , (3 Z = Ag x g y (53) 
The RG flows possess the following RG invariants satisfying Yli&dgiQ — 0: 

Qx = g 2 z - 9y, Q y = g 2 z - gl, Q z = g 2 x - g 2 y (54) 

There are only two independent invariants since Q z = Q x — Q y . We will thus express 
everything in terms of Q x , Q y . For the model defined by the action (jSOJ), the su(2) symmetry 
is maximally broken. However when g x = g y the U(l) symmetry generated by the currents 
J z , J z is preserved. This symmetry will guide us to classify the possible models. Observe 
that the self dual constraint q_i_ = Q-( Q i, = 0), discussed at the end of section III.B in 
connection with the models in |lRL \2(L |2l| , is not preserved by the RG flow unless g x = or 
g z = 0. If g x = the model is obviously gaussian, while if g z = the model is also gaussian 

n 

with exponents depending on the value of g x |2i|. These models where two couplings g's 
vanish simultaneously correspond to degenerate situations that shall not be consider in what 
follows. 

Let us study the RG evolution of the coupling g z associated to this U(l) symmetry. We 
can distinguish two cases: 



i 



gl>gl>gl Q x > Q y > (55) 



(ii) g 2 z <gl<g 2 y = : Qx < Q y < o 

The situation where g\ > g\ > g^ leads, in the U(l) limit, to the isotropic case and it is thus 
contained in the isotropic limit of (i) and (ii). Using the RG invariants one can eliminate 
g x ,g y from (3 Z : 



dg 



z 



(U As x s y ^{gl-Q x ){gl-Q y ) (56) 
14 



where (s x ,s y ) = (siga(g x ) , sign(g y )) . The solution of the above equation for the cases (i) 
and (ii) can be expressed in terms of Jacobi elliptic functions: 

9z(l) = \l~Q~x ns {^\[Q~xS x S y (loo - 1); &rg) , Qx > (57) 

9z(l) = V-Qx cs (y^-QxSxSyiloo - I); k' Ig j , Q x <0 (58) 

where is the value of I at which g z (l oo) — °°j an d ns(z; k) = l/sn(^;A;), cs(z; fc) = 
cn(z; ifc)/sa(z; fc) are Jacobi elliptic functions of z with modulus fcQ. ' The modulus of the 
elliptic function ns(z; fc) is here denoted as k Tg , and is defined as 

k 2 TS = < k% < 1 (59) 

while the modulus of cs(z; k) is the complement of k rg ( k' 2 g — 1 — fcj g ). The solutions ()57I58|) 
can be mapped into one another using the equation 

ns(ra; k) = —% cs(u; k'), k 2 + k' 2 = 1 (60) 

so that (|H7jl turns into (JH%|) by writing \[Q^ = zy '—Q x . 
The functions ns(z; fc), cs(;z; fc) have double periodicity 

ns(z + 4mK + 2mK'; fc) = ns(z; k) (61) 
cs(;z + 2mK + 4mK'; fc) = cs(z; k) 

where m,n are integers, K(&) is the complete elliptic integral, and K'(fc) = K(fc'), with 
fc' = a/1 — The coupling <^ is thus a periodic function of I, with a period depending on 
the sign of Q x : 

(i) Q x > 0. Here \[Q~ X is real and 

g,{l + \i) = -g,(l), Ai = ^^ (62) 

(ii) Q x < 0. Here a/QI is imaginary and 

g z (l + \[)=g z (l), K = ^7== (63) 

The 1-subscripts on A, A' refer to being the 1-loop result. The reason for the extra minus 
sign in eq. (|52*jl in comparison to eq. (J53J) will be explained below. 
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B. U{\) invariant limits 



In the U(l) invariant limit g 2 = g 2 , we have Q x = Q y = Q and k rg = 1. The solutions of 
the one-loop renormalization group again depend on the sign of Q: 



(i) Q > 0,k 



rg 



9z(l) = VQ coth [a^/Q - 0) (64) 

The model is well understood in this limit. Define \fQ to be positive. There are two 
subcases depending on the sign of g z . Massive case: When g z > 0, g z flows to an ultra- 
violet fixed point at g z = \[Q as I — > — oo and to a strong coupling fixed point in the 
infrared; it is thus a massive theory. Massless case: When g z < 0, the fixed point is at 
g z = —\[Q in the infrared (/ — > +oo). This is thus a massless theory. Both these theories 
have a sine-Gordon description where the conventional sine-Gordon coupling is related to 
Q; non-perturbative formulas were obtained in [(|, and will be used below. In this limit, 
since K(/c) ps log4/vl — k 2 as k — > 1, the period Ai in eq. (j62|l goes to oo, as expected for 
theories with fixed points. 

(ii) Q < 0, k rg = 1 Here since \J~Q is imaginary, as in Q| let us define: 



'Q = ih/A (65) 
The one-loop RG flow is now 

h 

g z (l) = - cot(/i(/ 00 - /)) (66) 

In this case the periodicity eq. (|63|) is maintained. Since K'(/c) pa 7r/2 as k — > 1, the period 
X[ becomes 

K - \ (k rg = 1) (67) 
This agrees with the manifest periodicity in eq. (J6TI|) and the one-loop period computed in 

Q 

In the U(l) invariant limit the higher order contributions to the beta functions are 

n 

known [24,]. These beta functions were used to classify the various phases as a function 
of \[Q in 0. The various regimes are distinguished by whether the fixed points are in the 
IR or UV, or whether the flow is cyclic. The result is summarized in figure 1. 
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FIG. 1: Different regimes in the U(l) invariant limit as a function of \[Q. 

In order to classify the possible models when k Tg ^ 1, we assume that each distinct regime 
of the U(l) invariant model has an elliptic deformation that is integrable. This is one of 
the main hypotheses of this paper. As we will show in the sequel, the fact that consistent 
S-matrices can be proposed for all these regimes supports this hypothesis. Our nomenclature 
for the models refers to the k Tg = 1 limit, and Q refers to Q x in this limit. We can thus 
identify the following distinct models: 

Elliptic sine-Gordon model (EsG). This is an elliptic deformation of the usual sine- 
Gordon model with < \fQ < oo. 

Elliptic cyclic sine-Gordon model (EcsG). This is an elliptic deformation of the cyclic 
sine-Gordon model described in [2| where Q < and y/Q = ih/4, with h > 0. 

Massless elliptic sine-Gordon model (mEsG ). Here —1/2 < -JQ < 0. This massless 
model is characterized by having an infrared fixed point in the U(l) invariant limit. 

Elliptic sinh-Gordon model (EshG). Here — oo < ^/Q < —1/2 and the model reduces 
to the usual sinh-Gordon model with one massive scalar particle in the U(l) invariant limit. 

It is important to note that under permutations of the labels x, y, z, the models can be 
mapped into each other, and this provides certain consistency checks of our results. Consider 
first the permutation: 

9x 9y — ^ Qx y k T gQ x , k rg > 1/^rg (68) 
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Under this transformation, the model is mapped into itself, which means that the period 
of the RG Ai, X[ should be invariant. This is readily checked using Re(K(l/k)) = kK.(k). 
Consider next the permutation: 

9 y <-> 9z => Qx -> ~Qx, k rg -> K g (69) 

This implies that the models EsG and EcsG are essentially the same at complementary 
elliptic moduli. In particular, when g 2 = g 2 , Q y = 0, and thus k Tg = is another U(l) 
invariant limit of EsG. Since the sign of Q x is flipped, this means that the EsG model at 
k Tg = can be mapped onto the EcsG at k rg = 1, and visa versa. The extra minus sign in 
eq. (I62|) was chosen so that the periods A, A' agree under this exchange. In the sequel this 
will serve as an important consistency check of the S-matrices. Though the two models EsG 
and EcsG are essentially the same, S-matrix descriptions given below are different because 
as defined the models have different U(l) limits as k rg = 1, one being the usual sine-Gordon 
model, the other the cyclic sine-Gordon model studied in |9(. However since the U(l) limit 
at k Ig = does not correspond to the conventional sine-Gordon action, but only after the 
permutation V(g z ) = g y ,V(g y ) = g z , as we will see the S-matrices in this limit match onto 
known U(l) invariant S-matrices up to a transformation V with V 2 = 1. Finally consider: 

9x ^ 9z — ^ Qx > k Tg Q x , k rg > i~j~ f (70) 

Here since k rg becomes imaginary, this does not lead to any equivalence between models. 

The above analysis is only at one loop and one must investigate whether the main features 
persist to higher orders. Our basic assumptions in the sequel are the following. First, the 
one-loop RG invariants Q, when appropriately corrected to higher orders, must continue to 
be RG invariants. Secondly, the RG flows must continue to be cyclic where the period of 
the RG, X(Q), is a function of the higher order corrected Q's. In the U(l) invariant limit, 



this is precisely the situation 



0,3 



When g 2 = g 2 one has: 



q - = q ' = a-f) 7 (i-a) i9 " = 3l) (71) 

The higher order beta-functions can also be integrated exactly in the U(l) limit, and the 
period of the RG in the cyclic regime is A = 7t/2a/— Q x , which is precisely twice the 1-loop 
result [dJ . 

In Appendix A, we study the higher order corrections and provide analytical and numer- 
ical evidence for the above hypotheses. There we give evidence that in the fully anisotropic 
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case the period is again twice the one loop result. This leads us to define the RG periods 
A, A' for the models EsG, EcsG respectively: 

A = 2Ai, A' = 2A; (72) 

where A^A^ as functions of Q x>y are given in eqns. (|62I63|) and the Q's are higher loop 
corrected ones. The fact that we can find exact S-matrices with the expected properties is 
a good indication of the validity of our hypotheses. 

In the above 1-loop RG analysis, the couplings are periodic on all length scales. The 
arguments given in section II suggest that all the elliptic models should be massless. This 
would however be inconsistent with some of the above U(l) limits which are known to be 
massive, unless a mass develops dynamically in going to this limit. It may also be that 
depending on what regularization one in practice uses to define the models, for example a 
lattice cut-off, the limit cycle may be only observable in the UV or IR. For this reason in the 
sequel we consider both the massive and massless cases. Most of the discussion will focus 
on the massive case since these results are straightforwardly extended to the massless case 
since Slr is the same function as in the massive case. (See section IX.) 

V. NON-DIAGONAL ELLIPTIC S-MATRICES 

In this section we review Zamolodchikov's elliptic S-matrix [jj], defering its relation to the 
quantum field theory to the next section. As usual, introduce relativistic massive dispersion 
relations eq. ([130. The S-matrices are more clearly presented in a real basis for the particles. 
Let A x , A 2 formally denote creation operators for the two particles in the theory. The S- 
matrix is encoded in the exchange relation: 

MPJMfo) = KtiPi - fh)Mfo)MPi) (73) 

The S-matrix £ does not possess a U(l) symmetry and thus has some additional non-zero 
amplitudes in comparison to e.g. the sine-Gordon S-matrix. Define: 
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(74) 
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Crossing symmetry in this basis reads: 

a((3) = a{m-(3), a t {f3)=a t (m-f3), a r {f3) = a a (m - /3) (75) 

The general solution to the constraints of the Yang-Baxter equation and crossing sym- 
metry has two free parameters, k, a: 

sn(2iria - 2a[3; k) 

<mp) = 777- — r\ — a( -P) 

sn(2ma; k) 
sn(2a/3;k) 

aM = sn(2^;fc) a{(3) (?6) 
cr t (P) = —k sn(2a/3; k) sn(27ua — 2a/3; k) a{f3) 

n 

(The parameter 77 in |2| is here expressed as rj = —ina.) 

When the S-matrix is real analytic, (S(/5))^ = £(—/?), unitarity reads E(/3)S(— /?) = 1. 
This gives the additional constraint on a: 

<r\ ( a\ sn 2 (27rm; k) 

a W°{-P) = m '(27ria;fc)-8n*(2a/?;Jfc) (7?) 



The so-called minimal solution to the above equation and crossing symmetry is 

1 sinh 2 (2irn(iT — sin (27m/3/7') sin {2im(m — f3)/j') 

n=l 

where 



j isinh (2nn(n - 7J/7 J sin (Zimp/j ) sin (27rn(wr - pj/7 J , > 

^—*n sinh (4^7/7') cosh (2ir 2 n/y) 



«W 2K<*) 
2q a 

The above S-matrix is real analytic so long as a is real. However the infinite sum in eq. 
()78|) is convergent only if 7 > tt/2 when 7' finite. For < k < 1, convergence of a thus 
requires that a be real and positive. 

In order to study U(l) invariant limits of the elliptic S-matrix, we go to a complex basis 
of particles: 

A± = (Ai ± zA 2 ) (80) 
The S-matrix for the A± particles is defined as 

MPi)Mfo) = S$(Pi - (3 2 - k, a) Mfo)MPi) (81) 
The non-zero amplitudes are 

C — C++ Q 

= — ^ 
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Q — CH — Q — + 

•J* — — *J_ _)_ 

S r = S++ = S+- (82) 

c — c — C++ 

°a — J ++ ~ J 



The relation between S and S follows from eq. (|50|l : 



Equivalent ly: 



2S = (7 — a + <7 t + ov 

25f = cr + a a + cr t - a r 

2S r = a + a a — a t + 0> (83) 

2S' a = cr — cr a — o-( — a r 



2a = S + S t + S r + S a 

2o t = S + S t - S r - S a (84) 



2a r — Sq — St + S r — S a 
2a a = —So + St + S r — S a 

Using (|75jl in eq.flHHJ) we obtain 

sn(27r«a — zap; re) 

*<ffl = =- ft(fl (85) 

sn(27rza — zap; re J 

5a(/3) = k sn(2a/3; re) sn(27ria; re) 5q(/3) 
where the modulus k and a are given in terms of the moduli k and a by, 

^(iT^) 2 - <H (1 + ^ a (86) 

For future reference, we note the identities: 

(1 + Vk) 2 K{k) = K'(re), (1 + V*fc) 2 K'(A;) = 4K(re) (87) 

The proof of eq.([85|) uses standard tools in the theory of doubly periodic meromorphic 
functions |4j. As an example let us consider the relation 

<j r Sq — St + S r — S a 



a S + S t + S r + S a 
21 



(88) 



which follows from eq.(fSlf). Using eqs.()76p and (jHHj) this equation becomes, 

sn(x — y; k) sn(x; k) — sn(y; k) + sn(x — y;k) — ksia(x; k)sia(y; fc)sn(x — y; k) 



(89) 



sn(x; k) sn(x; k) + sn(y; k) + sn(x — y; k) + /csn(x; fc)sn(?/; /c)sn(a; — /c) 
where 

x = 2nia, y = 2a/?, x = 2nia, y = 25/? (90) 

One can easily check that the LHS and the RHS of (J89|) . viewed as functions of x or y, 
have the same periodicity properties, position of poles and zeros, and hence by Liouville's 
theorem should be proportional up to a constant, whose value is actually one. 
In this basis, A + and A_ are charge conjugates and crossing symmetry reads: 



S t (P) = S (z7T - /?), SM = S r {m - P), S a (P) = S a {m - P) 



(91) 



It turns out that to describe all the field theories in section II, it is convenient to introduce 
a different description of the S-matrix. Define new particles A ± with the S-matrix exchange 
relation: 

AiWMfo) = Sf b {Pi - k, a) MP^MPJ (92) 
Define the non-zero amplitudes So,t,r,a as m ec L- (IH21), e.g. So = S+t = SZZ, etc, and let: 



So — <J r , St — cr a , S a — Ct, S r 



a 



(93) 



with the cr's the same as in eq. (|76p. Then it follows from the crossing symmetry relations 
on the cr's, eq. (j75]l . that the S"s satisfy the crossing relations eq. (|9~T|) with S — > S. Thus 
the S-matrix S is a proper S-matrix with A+, A_ charge conjugates. 
Expressing E as the matrix: 

(o a a \ 

(Jf (T r . 

o r o t 
\a a a / 

then in the sequel we will express S as the following transformation of S: 

/Sq S a \ /(T r <T t \ 



(94) 



S 



S t S r 
S r S t 
\S a So/ 



a a a 
(X a a 
\a t a r / 



(95) 
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where V 2 = 1. The V transformation can be written in matrix form as follows: 



S = P 2 Y>P 1 



(96) 



with 



where <r. 



Pi = Vol 
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1 
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(97) 



interchanges the two particles. Incidently, it is clear from eq. (|9T)|) 



that if £ satisfies the Yang-Baxter equation, then so does S. It should be kept in mind that 
the V transformation is not a change of basis. 



VI. MATCHING LAGRANGIAN AND S-MATRIX PARAMETERS FOR EsG 

As shown in section II, the field theories are symmetric and this should be a symmetry 
of the S-matrix. As shown in Q], this is indeed a symmetry of the elliptic S-matrices of the 
last section. The S-matrix parameters k, a are dimensionless parameters and thus must be 
RG invariant functions of the coupling constants g x ,y,z, equivalently functions of k Tg and Q x . 
By matching the periodicity of the RG with the periodicity of the S-matrix we now relate 
the S-matrix parameters k, a of the last section to the lagrangian parameters. In this section 
we do this for the EsG model. We assume the model to be massive; the massless version 
will be described in section IX. 



A. The S-matrix 

The S-matrices in eqs. (I73I74I93J) enjoy the following periodicity in rapidy: 

S(P-y)=S((3), S{fi-i) = S{p) (98) 

where 7' is defined in eq. ([79*)l. Based on the results of section III, we match this periodicity 
in rapidity with the periodicity of the RG and thus identify 

7' = 2A (99) 
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where A is the period of the RG defined in section II. The above equation implies 

K(Jfe) K(A; rg ) 



(100) 



Eq. ()100|) is a single equation for two parameters, however we can argue as follows to fix 
them both. When k rg is or 1, then the S-matrix must be trigonometric, i.e. k must also 
be or 1. However when k Tg = 1, k must also be 1, otherwise the relation between a and 
Q x would be infinite. This also implies that k rg = corresponds to k = 0. This suggests 
that k is only a function of k Tg , since both are between and 1. A further constraint is 
provided by requiring both U(l) invariant limits to be correct. Comparing the U(l) limits 
of the elliptic S-matrix to the usual sine-Gordon and cyclic sine-Gordon ones, this requires 
that when k = 1, a = y/Q/2 and when k = 0, a — \[Q. We have found the following 
solution to all these constraints: 

2\fk 

y/Q^=(l + k)a, h g = Y^ ( EsG ) ( 101 ) 
We have used the identity: 

K ) = (1 + k)K(k) (102) 

In summary we propose the S-matrix for EsG is: 

S EsG (P) = S(fck,a), (103) 

where S(/3] k, a) is given in eqs. (|82I83|) . and k, a are related to the lagrangian parameters 
Q x , k Ig by eq. (HOD). 



B. Resonance poles 

Above, we have related the S-matrix and lagrangian parameters by matching the peri- 
odicity properties of the S-matrix with the RG. In this subsection we show that the model 
also has the resonance poles with Russian doll scaling anticipated in section II. 

The function sn(z; k) has the following zeros and poles: 

sn(z; k) : zeros : z = 2n~K + 2imK.' 

poles : z = 2nK + (2m + l)zK' (104) 
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where m, n G Z- The zeros of sn(27rz5 — 25/5; k) thus lead to poles of S t>r at: 

(3 = — 7717/71") + 777' (105) 

where we have used the identity eq. (|87j). 

Requiring that there are no complex poles on the physical strip < Im(j3) < tt leads to 
the constraint: 

7 > 7T (106) 

Note that the latter is compatible with the convergence of a, which requires 7 > tt/2. When 
7 > 7T, there remain resonance poles on the strip — 7r < lm{0) < 0: 

(3 n = Vn- iVn, = 71"/, 7] n = 7 - TT (107) 

Here, since the model is assumed massive, as explained in section II, n must be positive. 
When 7 > 2ir there are no resonances since 7] n > 1. Finally, since 7' = 2A, one sees that the 
spectrum of resonances precisely satisfies the Russian doll scaling property eq. ([22)1 . That 
the S-matrix has both of the UV signatures of a cyclic RG described in section II, with 
compatible period, relies on the special relation between the periods and poles enjoyed by 
Jacobi elliptic functions. 



C. U(l) invariant limits 

Checks of the above S-matrix are the U(l) invariant limits k vg = 0, 1. This check is non- 
trivial because the lagrangian and S-matrix parameters were related by only matching the 
periodicity of the S-matrix, and in the k rg = 1 sine-Gordon limit this periodicity disappears. 
The limit k rg = on the other hand preserves the cyclicity of the RG. As we now show the 
S-matrix gives the expected limit in both cases despite the fact that these two limits have 
very different properties. 

We first begin with k rg = 1, which is expected to be the sine-Gordon model. The latter 
is defined by the lagrangian 

S = -L J d 2 x Q(<90) 2 + A cos(60)^) (108) 

with < b 2 < 2. In the limit k Tg = 1, Q x — Qy = Q, and k = 1. The relation between 
b and Q can be found by matching the slope of the beta function at the fixed point and is 
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known [ 

2 V^ = 2 -y^ (109) 

Note that in this limit the sine-Gordon coupling g depends on both g z and g x , not only on 
g z as one would naively think. The reason is that both g z and g x continue to flow under 
RG. On the other hand b, being a dimensionless parameter in the S-matrix, must be an RG 
invariant, as in eq. (jlU9|) . This was explained in more detail in fl. In , 7' — > oo, 

and thus the RG period A becomes infinite, consistent with a theory with a fixed point. The 
resonances become infinitely heavy and decouple. 

The parameter 7 remains finite and depends on the sine-Gordon coupling constant b: 

71 71 7x1)2 mm 
7 = 4^ = 27^ = ^ (U0) 

Since 7' = 00, the constraint 7 > tt is not required and the whole range of the massive 

sine- Gordon model < b 2 < 2 is covered. 

Since k = in the limit, using sn(z; 0) = sin z, one finds from (|53|l S a = and 

s = sinh(4/fa) s 
sinh(4a(i7r — (3)) 
sin (47m) 

S r = i . rr A y. ^tt S (111) 

smh(4a;(z7r — p)) 

sinh(2a(i7r — (3)) 

cosh(2a/5) sinh(27ria) 

Since 7' goes to 00, the expression for a leads to an integral. Using in addition the integral: 

f°° dxcosh(ax) — 1 

logcos(2aa) = - / . , ; , (112) 

J x sinn(7rx/4a) 

which is valid for |o| < vr /4a, one can represent the additional trigonometric factors in eq. 
fllllj) with c = it /4a, and one finds: 

J x cosh(7rx/2) smh(7ra;/8a) 
Using eq. (|110|) . this agrees with the known sine-Gordon S-matrix[l|, with precisely the 
right dependence on the sine-Gordon coupling b. (So can also be expressed as an infinite 
product of T functions; see appendix B.) 

Consider next the limit k rg = k = 0, with Q x = Q. In this limit, 

*=VQ = j (H4) 
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and the period A remains finite: A = lnjh. This theory should correspond to the cyclic 
sine-Gordon model of jjj. In order to compare this limit of SesG with the S-matrix in 3], it 
is necessary to go to the A\^i basis of particles with S-matrix E. When k — 0, the non-zero 
amplitudes are a, a r , a a . As explained in section IV, in order to compare with the usual cyclic 
sine-Gordon model, which is defined when g x = g y , one must make the V transformation 
which exchanges g y and g z . Making this P-transformation to S defined by eq. (|95|h one 
finds in this limit: 

? sin(/^/2) ? 

* sin(/i(i7r-/5)/2) 
- , sinh(7r/t/2) g 
sm(/i(z7T — p)/2) 

S ° ~ ~ % sinh(vrV2) fffc=0 
In this case, since 7 — ► 00 and 7' remains finite, the expression eq. ()78|) becomes an infinite 
sum rather than an integral: 



n=l 



- zlog£ (/3) = vr + W2 + £ I 7 (n f } .. (116) 
The above S-matrix agrees with the one in 0|. 



VII. S-MATRIX FOR EcsG 

As explained in section IV, the EsG and EcsG models are essentially different descrip- 
tions of the same theory at complimentary elliptic moduli. The discussion therefore closely 
parallels that of the last section, and so we provide less details. 

A. S-matrix 

As we now argue, the S-matrix for this theory is more naturally described by S. As for 
EsG, matching the period 7' of the S-matrix with the period of the RG A', 7' = 2A', leads 
to: 

K(jfe) _ iK'(fc rg ) 



117) 



Kg = TTT V = &=(l+k)a (EcsG) (118) 



Repeating the arguments that led to eq. (jl()l|) one obtains 

1 + k 
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Note that since Q x < 0, a is still real. In summary the S-matrix is 

S BcsG (f3) = S(P;k J a) J (119) 

where S(/3; k, a) is defined in eqs. f)93jl . and the parameters k, a given in terms of Q x , k vg in 
eq. (ITT51) . 

B. Resonance poles 

The poles of sn(2iria — 2af3) lead to poles at (3 = — (2m + 1)7/71") + wy' /2 in the 
amplitudes Sq . As for EsG, requiring that there are no complex poles on the physical strip 
leads to the constraint 7 > tt. Given that the above constraint is satisfied, there remains 
resonance poles at: 

Pn = Hn-iVn, /J>n = «t7 2 > Vn = 7 - TT (120) 

Again, for the massive case n must be positive and these resonances only exist for 7 < 27r. 
Noting that — /i„ = 7'/2 = A' the Russian doll scaling of the resonances in the UV 

is: 

M n+2 w e v M n , for n large (121) 

This is still consistent with the UV signature of a cyclic RG described in section II, eq. (|18jl. 
the difference being that two states are reshuffled in each RG cycle A'. 

C. U(l) invariant limits 

Again the U(l) invariant limits k rg = 0, 1 serve as non-trivial checks. First consider 
k Tg = 1. In this limit, k = 0, and we parameterize \[Ql = \fQ = z/i/4 as in eq. (|65|1. One 
finds in this limit 

a = = h/A (122) 

and the period of the RG remains finite, A' = 2ir/h, and agrees with the RG calculation in 
section II. 

The parameter 7 on the other hand becomes infinite and the resonances disappear from 
the spectrum since 77 > 1 and they are not on the strip — tt < Jm(/3) < 0. Note also that 
since K'(0) = 00, the constraint eq. (jl(J6J) is not required and < h < 00. 
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One finds for the S-matrix the same result as in eq. ()115|) , and again this S-matrix agrees 
with the cyclic sine-Gordon one in 

Finally we consider the k Tg = limit. Here the period of the RG goes to infinity and the 
theory should be equivalent to the sine-Gordon model with coupling b and Q = —Q x related 
as in eq. (jl(J9|) . From eq. (|118|) . this limit corresponds to k = 1 which leads to a = y/Q/2. 
In order to make contact with the usual sine-Gordon S-matrix in this limit, we first make 
the transformation to X = V(S), then make a change of particle basis so that the S-matrix 
is given by S in equations (|83|85|) . Using now k = 0,5 = 2ia, one finds the result in eq. 
(jlllj) . with again correct dependence on the sine-Gordon coupling b. 



VIII. SCALAR THEORIES: ELLIPTIC SINH-GORDON 



The sinh-Gordon model is defined by the action 

S shG = ^J d2x (\( d( t>? + A cosh bcj^j (123) 

As explained in this model is realized in the field theory of section II when Q x = Q y = Q 
with ^fQ < —1/2. The relation between Q and b is known j^]: 



2VQ = -{2 + b 2 )/b\ (124) 
which follows from eq. (jl(J9|) with b — > ib. 



The spectrum of the model consists of a single massive scalar particle with S-matrix 



tanh(/3 - i7ra)/2 
tanh(p + tira) / 2 

where 

a = ^ («) 

Mussardo and Penati have considered the simplest possible scalar S-matrix built out of 
elliptic functions j^J: 

q(R u n \ - sn(2K^/7r;fc) + sn(2Ka;A;) 
lP ' ' >~ sn(2Ki/3/7r; fc) - sn(2Ka; fc) 1 ' 

where K = K(/c). In the limit where the elliptic modulus k — > 0, one recovers the sinh- 
Gordon S-matrix eq. (|125j) . 

We propose that the above elliptic S-matrix describes the elliptic sinh-Gordon regime 
of the field theory in section II. As before, we relate the S-matrix parameters k, a to the 
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lagrangian parameters Q x ,k Tg by matching the periodicities. The above S-matrix has the 
periodicity: 

S(p-2\;k,a) = S(J3;k,a), A = - ^ (128) 

Identifying A as the period of the RG, which from section IV equals 2Ai since Q x is positive, 
one obtains: 

K(k rg ) _ nK'(k) 



VQ~ X 2 K(k) 

Since the equation (j!29|) does not depend on a, we fix a using eqns. (jl24H26|) : 

' (130) 



2 yQx 

The above identification guarantees that the S-matrix has the correct limit as k — > 0. Given 
the above identification of a, then the value of k is determined by eq. 



For < A; rg < 1, a > 0, there is always a solution of the above equation with < k < 1. 
This completes the identification of k, a in terms of Q x , k rg . Near A; rg = 1 one has: 

k w 4(*4/4) 2a (132) 

so that approaches as fc rg — > 1. 

The above S-matrix has no complex poles on the physical strip as long as < a < 1. In 
terms of the lagrangian parameter this reads 

y/Q~ x <-l/2 (133) 

In the usual sinh-Gordon limit k — > 0, Q x — > Q, the above constraint correctly goes over 
to the sinh-Gordon regime (see section II), which provides a check on the S-matrix and eq. 

(unnD- 

The remaining resonance poles are at: 

Pn = fJ> n - i>Vn, V-n = 2nA, r) n = ira (134) 
Again these resonance poles satisfy the expected UV Russian doll scaling eq. 
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IX. MASSLESS ELLIPTIC SINE-GORDON MODEL (mEsG) 



The only region not covered in the U(l) invariant limit by previous cases is —1/2 < y/Q < 
0. The U(l) invariant model has an infrared fixed point and is thus massless. When the 
model is fully isotropic, y[Q = 0, it corresponds to the large distance limit of the 0(3) sigma 
model at 9 = 7rj3jJ. In this section we propose an S- matrix when k Tg ^ 1. The discussion 
closely parallel's the EsG case of section V. 

An S-matrix description of massless theories was given by Zamolodchikov and 



Zamolodchikov 31 1. An essential ingredient of their formulation are S-matrices for only 
left-movers or only right-movers, denoted Sll and Srr. These S-matrices are formal S- 
matrices for a scale-invariant theory. When left-right scattering Sm, is non-trivial, the scale 
invariance is broken. 

Since \fQx~ is rea l> t ne period of the RG is A = 2X%. Matching the periodicity 7' of the 



S-matrix while requiring a > gives: 

(mEsG) (135) 



K(k) _ K(fc rg ) 



The RHS is positive since here \[Q~ X is negative. As for EsG, we argue that the above 
relation requires: 

y/o;=-(i + k)a, k ^ = Y^k (136) 

Let us parameterize the massless energy momentum for left and right movers as in eq. 
(|24|). Requiring the two-particle S-matrices to correspond to the 0(3) sigma model at 9 = tt 
in the su(2) invariant limit fc rg = 1, Q — leads to the obvious proposal: 

Srr(J3) = S EsG ((3), P = Pm-Pm 

SlM = S EsG ((3), P = Pli- Pl2 (137) 

SrM = S EsG (P), P = P m -PL2 

where all the S-matrices 5*e s g o n the right hand side are the same as in eq. , where now 
a, k are determined by eq. (|136|) . As for other cases, one can easily check that this has the 
correct limit as fc rg = 1. 

The S-matrix has the periodicity eq. (J26j) . The abscence of complex poles on the physical 
strip requires 7 > tt. Using eq. (|136|) this gives \fQ~ x > —(1 + k)K'(k)/2n. Since K' > x/2, 
near k Tg = 1, this gives \[Q > —1/2, which is the expected range. 
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The analysis of resonance poles closely parallels the discussion in section VI. The poles 
eq. (|107|) are still present, but now as explained in section II, /i n can be negative. This leads 
to a spectrum of resonances: 

M n = m e nX cos((7 - tt)/2) n G Z (138) 

where the above equation is exact since the theory is massless. The above resonances 
accumulate at zero as n — > — oo and are infinitely heavy as n — > +00, which are both the 
UV and IR Russian doll signatures. This theory is thus consistent with a cyclic RG on all 
scales. 



X. CONCLUSIONS 



In summary, based on the limit-cycle behavior of the RG for maximally anisotropic 
su(2) current interactions, we have proposed that they underly the previously known exact 
S-matrices built out of elliptic functions. The S-matrix and lagrangian parameters were 
related by matching the period of the RG with the UV signature of a cyclic RG in the 
S-matrix, i.e. the periodicity in rapidity. Numerous checks were performed, in particular 
we showed the models have an infinite spectrum of resonances consistent with the cylic RG 
and have shown that the models have the expected U(l) invariant limits. 

Since in this paper we have proposed a field theory for the elliptic S-matrices for the first 
time, there are many open avenues for further investigation, and we finish this paper by 
listing a few of them. 

The sine-Gordon theory, which appears in the trigonometric limit of our model, is clas- 
sically integrable, and in fact semi-classical methods using this integrability were used early 
on|32| to study the spectrum of the model and these results eventually provided some checks 
on the S-matrixil]. To our knowledge the classical (and quantum) integrability of our field 
theory has not been studied. Clearly we have assumed it was integrable in proposing exact 
S-matrices. A good starting point is the bosonized action (|38jl . 

The sine-Gordon theory which arises in the trigonometric limit is known to have a quan- 
tum affine sl(2) q symmetry and the conserved charges can be constructed explicitly in the 
quantum field theory 3^ . It would be interesting to extend this field theory construction 
to the present models since this should lead to an elliptic deformation of the affine sl(2) 
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algebra which can be compared with the algebra in 3^ . 

Our model may have some applications to solid-state physics. In the fermionic represen- 
tation of the su(2) currents, the model is essentially a Luttinger liquid for fermions with 
spin and 2 additional kinds of density-density perturbations since the usual Luttinger liquid 
corresponds to g x = g y = 0. 

Theories with an infinite number of resonances are reminiscent of string theories j^]. In 
string theory the resonances are exactly stable, whereas in our model they generally have 
a finite lifetime. Another S-matrix was studied in which, though related, is essentially 
different from the S-matrices here, and is characterized by an infinite number of exactly 
stable resonances but with no periodicity in rapidity. A field theory interpretation of this 
S-matrix seems unlikely since it suffers from a lack of real analyticity. How this S-matrix is 
related to the physical S-matrices in this paper is described in appendix B. 
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XI. APPENDIX A: HIGHER ORDERS 



In 



24j an all orders beta function was proposed for general anisotropic current interac- 



tions. In addition to the coefficients C defined in eq. (JHOj) the beta function depends on the 
following. Expressing the perturbing operators as O a = d A b J a J b , define the purely chiral 
operators 

T\z) = diJ\z)j\z) (139) 
In the conformal field theory one has the closed operator product expansion 



T A (z)O B (0) ~ -V [2kD£ s + C£ B ) O c (Q) 



(140) 



The formula in 



24j for the beta function is then expressed in terms of the coefficients C,C,D. 
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For our model one easily finds the non-zero values: 

D xx = D yy =D zz = 2 (141) 
gxy = (jxz = (jyx = jjyz = (jzx = jjzy = 4 

From the formula given in [2^ one finds 

= ^9x9y(l + k 2 9 2 z )-kg z (9 2 x + 9y)) (U3) 
(l-K»<g)(l-£flg) 

where here k is the level of the current algebra, which for our model equals 1. The other two 
beta functions /3 X , j3 y follow from the above expression by permutation of the x, y, z indices. 

The above beta function has a strong-weak coupling duality. For the U(l) invariant case 
9x = 9y-> this duality was exploited in [6[ in order to extend the flows to all scales. Since the 
level can be easily scaled out of the equations, let us set k = 1. For each coupling, define 
the dual coupling as g* = 1/g, and the beta function for the dual couplings: 

W)=i3(j)y (144) 
Then one can verify that the beta function satisfies the duality relation: 

P*(g*) = -P(g^9*) (145) 

An important consequence of the duality is that flows between g equal to and 1 can be 
mapped onto flows between g equal to 1 and 00. We can now use this duality to argue that 
the flows based on the above all-orders beta functions continue to be cyclic as follows. At the 
intial RG time, suppose that the couplings are near zero, and that running forward in RG 
time they reach g = 1 whereas running backwards they reach g = — 1. Though the points 
g = ±1 are poles in the beta functions, it was shown in the U(l) invariant case that because 
of the RG invariants the flows approach the poles along tangent directions determined by Q 
and flow smoothly through the pole. In other words, the poles are not true singular points: 
a local blow up resolves the flows. Beyond g — 1, the flow between g = 1 and 00 is dual to 
the flow between and 1. Thus it takes the same time to flow between 1 and 00 as it does 
between and 1. At g = 00 the flow actually continues smoothly at g = —00; this jump is 
exactly dual to a smooth flow through a g = since I/O 1 * 1 = ±00. The flow then continues 
to g — — 1 and a new cycle begins. The period of the RG is then twice the time it takes to 
flow between —1 and 1. We will use this below to numerically determine the RG period. 
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What complicates the situation in the fully anisotropic case is that unlike the U(l) invari- 
ant one we have been unable to find simple expressions for the higher order RG invariants 
Q, nor have we been able to integrate the all orders beta function and compute the period 
A analytically. However we can at least show that the existence of the RG invariants is not 
spoiled at 2-loops. Additionally, the two-loop analysis shows that the corrections to the Q's 
are not simple power series. Let (3(g) = where (3^ is the i-th loop contribution 

the beta function, of order g l+1 . Let us also write for Q = Q x , y ,z' 

Q = Q« + + . . . . (146) 

where Q^ 1 ' is the one-loop expression given in eq. (J54j) . and Q is the two loop correction. 
Assuming is of order g 3 , the RG invariance of Q to two loops requires 

£ pPd gi QW+(3Pd 9i QM=0 (147) 

i=x,y,z 

Keeping just the two-loop contributions to eq. (|143|) one finds 

(3 z = Ag x g y -Ag z (g 2 x + g 2 y ) + .... (148) 

and (3 x>y again given by the obvious permutations. Using this in eq. ()147J) one finds the 
following two loop correction to Q x : 



Qx = (g 2 z ~ 9 



2\ 



1 + 2 



ll.r. 



u 2 du 



y V + g 2 y - gl){u 2 + 



9z ~ 9 2 x ) 



+ .... (149) 



The above is an elliptic integral of the third kind. This suggests that the RG flows can be 
uniformized using elliptic functions, reminiscent of Seiberg-Witten theory [sfij]. 

We now can give numerical evidence for one of the main hypotheses of this paper, i.e. 
that the period of the RG is twice the 1-loop result, eq. (|T2*j) . If the couplings are initially 
very small, then we expect that we can approximate the Q's by their one loop expressions 
Q^K The results in the case Q x < for several values of initial couplings are shown in 
Table 1. The analytic expression A' = 2A' 1; is given by the formula (|63|) . where Q x ,Q y are 
approximated by the 1-loop result Qx}y The last column, A num , denotes the period of the 
exact RG evolution using the beta functions eq. (j!43jl . which we have computed numerically 
as the time to go from a g z = — 1 to g z = 1. In the U(l) case, where g x = g y = 0.1, there 
is a small discrepancy between 2X[ and A num which is an indication that we have used the 
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1-loop expressions in the expression for Ai. This presumably explains also explains the 
small discrepancy when k rg ^ 1. 



9x 


9v 


9z 


h 


Qi 1} 


Qy 


2A' X 


Anum 


0.1 


0.1 





1 


-0.01 


-0.01 


15.708 


15.629 


0.1 


0.2 





0.5 


-0.04 


-0.01 


10.783 


10.660 


0.1 


0.3 





0.333 


-0.09 


-0.01 


8.429 


8.2590 



Table 1.- Numerical comparison between the analytic expression 2X[ and the numerically 

determined period A num . 



XII. APPENDIX B: INFINITE PRODUCTS OF T-FUNCTIONS AND A 
STRINGY S-MATRIX 

The other S-matrix considered in [^] was an analytic extension of the usual sine-Gordon 
one to the complex values of the sine-Gordon coupling b: 

/— 2 - b 2 ih 

n 

As pointed out in |9|, the resulting S-matrix is not real analytic: S^((3) ^ S(—/3), and for this 
reason more than likely does not have a field theory description. Regardless, the S-matrix 
has some interesting properties. It is characterized by a Russian doll spectrum of resonances 
in the UV but with no periodicity in rapidity. The resonances are exactly stable, and closing 
the bootstrap led to a string-like spectrum 9J]. 

This stringy S-matrix has not played any role in this paper, nevertheless it is closely 
related to the EsG S-matrix with a purely imaginary, as we now explain. Consider the 
EsG S-matrix in the limit k — > 1, which leads to eq. (jlllj) . Let a be defined by eqns. 
(111011501) : 

a=- (151) 
8 V ' 

Then the S-matrix in eq. ()111|) has the overall structure of ratios of the cyclic sine-Gordon 
model in eq. ()115j) with S — > S. When k — > 1, 7' becomes infinite and 7 remains finite, thus 
the expression for a eq. ()78j) becomes an integral. However the integral does not converge 
for a purely imaginary. 
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In |9j a different overall scalar factor was proposed which is just the analytic extension 
of the usual sine-Gordon one. Define: 

cr _ r(4a)r(l - z) -ft F n (0)F n {iir - g) 

6 ° " T(4a-z) 11 F n (0)F n ^) (152) 

where 

r (8na - Z )T(1 + 8na - z) 
n{ P) r (4(2n + l)a - z)T(l+ 4(2n -l)a-z)" 1 J 

and we have defined: 

4i/3a 



(154) 

7T 

When a is real and related to the sine-Gordon coupling b as in eq. (jllUj) . the above is 
the well-known expression for So as an infinite product of T-functions. (We will show it is 
equivalent to eq. (|113|) when a is real below.) 

Consider now the above infinite product when a is pure imaginary. The product is still 
convergent. It can be given an integral representation using: 

Ite'" 81 -" ' l-expt-^/to) = l0gr(4a0) + ^ ~ ll2){ ~< ~ ' 0g4a) " l0g(2T)/2 

(155) 

where 7 is Euler's constant, and the contour C is shown in figure 2. The above integral is 
valid for a real or purely imaginary as long as the real part of a positive. Nearly all the V 
functions in logS^ can be represented with the real part of a a positive integer. The sum 
over n converges since Y2n>o ex P( — 2n7rx) converges. The result is: 

, r T(l + Ai(3a/7i)T(Aa-Aif3a/7T) 



where 



r(l-4z/3a/7r)r(4a + 4^a/7r) 1 ; 

log/- / ^ lQg( Tlx) sin (^) silin ( 7rX ( 1 ~ V 4 ^)/ 2 ) c -ttx ^ 57 ^ 

J c 2ttx cosh(7rx/2) sinh(7rx/4a) 

In eq. (|156jh we have factored out the T functions that cannot be represented by an 
integral when a is imaginary. The above integral is convergent which proves the infinite 
product is convergent. 

When a is real, the additional F functions in eq. (|113j) can also be represented by a 
contour integral. Furthermore, since in this case there are no poles on the real x axis, the 
contour integral can be replaced by an ordinary integral: 

/ -^log(-vrz) — (158) 
Jc 2nix J x 
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|_x 



FIG. 2: Integration contour C for eq. ()155|) . 



and one recovers eq. (|113j) . For a imaginary however, because of the 1/ sinh(7ra;/8a) in the 
integrand one cannot safely make the replacement eq. (jl58j) . 

In summary, when a is imaginary, one cannot obtain the convergent expression eq. (|156|) 
from the k — *■ 1 limit of a in eq. (|78j) . Rather, one has to perform the integrals in eq. piHjl 
for a real obtaining the infinite product of T functions, and then analtyically continue to 
imaginary a. 
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